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^— ■») ' We investigate the effect of an incipient ferrolectric transition on vertex corrections to the super- 

\ conducting pairing interaction. The vertex corrections for small momentum transfers are large inde- 

, ^ . pendent of the type of Boson responsible for the superconducting transition. The electron-phonon 

^ interaction is found to be enhanced by a nearly ferroelectric medium. We discuss application to the 

, cuprate superconductors. 

in 

(N 

1) INTRODUCTION. 

^ ■ High-temperature superconductivity is commonly believed to be associated with the exchange of antiferromagnetic 
O spin fluctuations. A unified description of antiferromagnetism and superconductivity is given by a S0(5) group scenario 
^ 0. However, it had been previously shown by Birman and Solomon 1^ that a similar group-theoretical formulation is 
^ possible with superconductivity associated with charge-density waves, ferromagnetism, and ferroelectricity, as well as 
the "common" antiferromagnetic scenario. We consider here an effect associated with the nearness of a ferroelectric 
^ ' transition. A phenomenological pairing interaction has been proposed for such a nearly ferroelectric system and 
^ ' one purpose of this paper is to make a first step towards giving this pairing interaction a firmer microscopic basis. 

Quite generally, the pair interaction can be viewed as the exchange of a Boson between the members of a Cooper pair. 
^ Of primary interest here is the electron-Boson vertex. We re-evaluate the Bardeen calculation of the correction to 
the electron-phonon vertex due to electron-electron interactions, for the case of an electron-gas imbedded in a dielectric 
medium with a large dynamic dielectric constant e reflecting the nearness to a ferroelectric transition. We calculate 
the lowest order vertex correction that includes phenomenologically the important ionic dielectric effects. This turns 
^ out to be a large correction. Although we assume in the following that the Boson is a phonon, it is important to note 
that the present result is not limited to this case and indeed would also yield a large vertex correction in the popular 
scenario in which the Boson is an antiferromagnetic spin fluctuation. 

Many perovskites arc fcrroelectrics. Perovskites like SrTiOa and the high-Tc cuprates are nearly ferroelectric; their 
ionic dielectric constant is abnormally high, although they do not undergo a ferroelectric transition. The dielectric 
' constant of La2-a;Sr2;Cu04 and YBa2Cu307_5 was measured recently quite accurately as function of frequency and 
. is of the order of 100 at frequencies of order 10 meV |^J^. In "conventional" calculations of the electron-phonon 
' interaction, this dielectric constant is not taken into account. The reason is that, although the electron-phonon 
, interaction itself is a nonadiabatic contribution, it is calculated within the framework of the Born-Oppenheimer 
1^ ■ approximation (BOA); i.e., it is assumed that during the scattering of the electron, while it imparts momentum to the 
' lattice, the lattice does not move. In contrast, the ionic dielectric constant describes motion of the lattice. Therefore, 
I in order to take it into account, we must consider effects outside the BOA. This presents a fundamental difficulty, 
' and not just "technical" considerations, 
i' . The phonon-mcdiated interaction considered by Bardeen and Pines in the 1950's involves an effect outside the BOA. 
'■^ ' The scattering of electron 1 causes an atomic motion, and this motion interacts with electron 2. The Bardeen-Pines 
^ i interaction is given by : V = g^ ui^-n^ ' where f2 is the phonon frequency, g is the Frohlich constant given by g = 
Q ' I/V MVl, M is the ionic mass, I is the matrix element given by I(q) — {k \ \ k + q) , and V is the electron-ion 
• • potential. I(q) is thus calculated within the BOA. This is an apparent paradox; this procedure is justified by Migdal's 
theorem which states that the corrections to I(q), arising from the fact that we consider an effect outside the BOA, 



' are of order fl/Ep ~ yJmjM ~ 10^^ , where m is the electronic band mass. 
^ ' Pietronero et al ^ suggest that in the high-Tc cuprates fl/Ep is not so small, and as a result the deviations 
from Migdal's theorem are significant. Here, we also question the validity of the BOA, but for a different reason - 
namely, we suggest that the very large ionic dielectric constant has a very large effect on the electron-phonon matrix 
element I(q). It renormalizes it, and thus makes it frequency dependent; at very low frequencies I(q) is increased by 
a significant amount above the BOA value. 

2) CONVENTIONAL APPROACH. 

The "conventional" way to calculate the electron-phonon matrix element I(q) is based, first, on Bloch's calculation 
of 1928 0, in which he assumed: I(q) = (fc | WV \ k + q) , where V is the electron-ion potential Ze^/r . Thus I(q) 
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= AwZe^/q . I(q) is seen to diverge at small q values. We denote this as the "Bloch Vertex" TBioch (Fig. la). This 
estimate neglects the electron-electron interaction which screens the electron-ion potential. This screening was first 
considered by Bardeen in 1937 [||. The electron-electron interaction propagator is given by: Dee (q) — 47re^/g^ . 
Bardeen took it into account by considering the potential: 



Vsardeeniq) - i , n / „\ AT ^ \ ^ ^2 i „2 ' i^) 



q^ — ATTe'^N{Ep). This potential no longer diverges as q— > 0. The electron-phonon matrix element is given by: 

I(c[) = ^^2^^]'' • Since qu > , the dependence of I(q) on q is weak and the electron-phonon scattering is nearly 

isotropic [0. The Cooper potential reflects this isotropic nature and leads to superconducting pairing with s-wave 
symmetry. What Bardeen did to deal with this (without saying it explicitely) was to sum-up RPA bubble diagrams 
of the form of Fig lb to get the correct screened vertex. We call this the "Bardeen vertex" V Bardeen ■ 

The simplest vertex correction arising from the electron-electron interaction, /^Tcoulomb^ is shown in Fig. Ic. We 
calculate this contribution in the next section for comparison with the vertex correction including ionic effects. In 
Fig. Ic we could of course replace Dee with a screened interaction of the form of Fig. lb. Note that Migdal's theorem 
does not apply here since Dee does not possess a low-energy cutoff. 

3) VERTEX CORRECTION INCLUDING THE IONIC DIELECTRIC FUNCTION 

In a medium with a dielectric constant e, the electron-electron interaction is given by: Ane^ /q^e . Here e is the 
dielectric constant of the medium external to the electron gas pT[| . Thus, we use a dressed electron-electron propagator: 

{q'^ + qfj)e{q,uj) 

Diagramatically, we denote Dee by a heavy broken line (Fig. Id) that denotes the inclusion of the effect of the highly 
polarizable ions on the electron-electron interaction. The Thomas Fermi screening vector qo is included to facilitate 
the numerical computations and to account for electronic screening. Thus we consider the vertex correction AF 
shown in Fig.le. For simplicity in the present calculation we take TBardeen to be a constant. We can thus write F = 
Fsardeenil + ( AF/Fsardeen)] and wc Calculate AT /T Bardeen and AT Coulomb Bardeen which are to be compared with 
the number 1. 

We have to calculate the vertex correction as a function of q, paying particular attention to small q values. This is 
plausible for several reasons. First, the ionic dielectric corrections could also be included in the Bardeen screening by 
replacing Dee with Dee and setting qjj = to avoid double-counting of screening diagrams. Then if Dee were zero (as 
is the case for an infinite e), F would be the Bloch vertex which diverges as q^ 0. Thus AT /T Bardeen would diverge 
as q^ 0. 

The necessity to consider small q values was also pointed out by Pietronero et al for a different reason. They 
show that for their (phonon) vertex correction, the behavior when q—^ first, and w — > afterwards, is entirely 
different from the behavior when oj ^ Q first, and q^ afterwards, which was previously calculated by Grabowski 
and Sham |l|. 

The relevance of small q values for the cuprates has also been pointed out previously by Santi et al |l3|, P erali and 
Varelogiannis [Q, Abrikosov ||l^, Zeyher and Kulic jl^, Bulut and Scalapino ||l^, Tsuei and Kirtley |18|, Leggett 



and others. The observation of stripes |20| suggests the relevance of q values close to q stripe — 0.24 A ^ 



and we calculate here the vertex correction for q values in this range. 

The dielectric constant has a sharp dispersion at a frequency of about 19 meV in YBCO Q. This frequency is 
the frequency of phonons involving the displacements of Barium atoms in the c-dircction. In LSCO , the dispersion 
frequency is about 27 meV, the frequency of phonons involving motion of strontium in the c-direction [Q. A theoretical 
expression for the ionic dielectric constant is given by the Lyddane-Sachs- Teller theory (LST) |2^, namely: 

2 , ,2 



UJ ~ OJt 
— U!^ 

where ul and ut are the frequencies of the longitudinal and transverse phonon modes. Although this expression is 
derived for q=0, we assume that it is valid for the small values of q of interest here. 
The use of such a simplistic expression, however, raises several further questions: 

(a) There are several phonon modes, and there is no a— priori reason why e{uj) should be dominated by one 
transverse, and one specific longitudinal mode. 
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(b) We ignore the dispersion of the phonon modes for the relatively small values of q that we consider here (about 
1/4-1/3 of the way to the Brillouin zone). This is justified here since there is no softening of the modes, as manifested 
by absence of a strong temperature dependence of their frequency. |^ 

(c) This expression applies to insulators. The cuprates possess a metallic Fermi surface. 

(d) The expression is derived for an isotropic 3D system. 

We feel reasonably confident using Eq.(3) simply because experiment shows that it wor ks quite well. §,§. For 
example, from optical infrared measurements, it is found that the simple LST expression applies extremely well to 
the c-axis component of e for LaSrCuO (with =63 meV ) 0, and reasonably well for YBCO [||; (there are two 
modes that serve as lol , at 40 meV and 70 meV; we can choose some average in-between these two frequencies; about 
50 meV gives good agreement, with — 4). 

Since the "metallic" layer is thin, the width being about 2ac where a.c is the Bohr radius of the oxygen 2p(T orbitals 
in the c-direction, &c ~ 0.4 A, and the average distance between electrons is approximately the lattice constant (~ 4 
A) , the c-axis component of e screens out the electron-electron interaction effectively, and we can use it (as determined 
experimentally) in the expression for D^e- 

The contribution AF, Fig. le, is given at zero temperature by 

AT{k,q)=TBardeenY. j ^ D,,{p- k)Go{p)Go{p + g) , (4) 

where k = (k,fco) and we assume a two dimensional system. Substituting Eq. (3) for e(a;) with eoo = 1 in Eq.(2) we 
have 

0,2^^1) kl-u;l + ^5■ 

This form is quite plausible but, as mentioned above, has not yet been derived from a truly microscopic theory. Go 
is the usual zero order electron propagator and we assume for now a free particle spectrum. The po integration can 
be caried out as a contour integration picking up contributions from the poles of the Go's and D^e- 

Since F is normally part of a larger diagram, the pair interaction, for example, its external momenta and energies are 
integrated over. Although for such applications we require the full complex AF as a function of all of its arguments, 
for present purposes it should suffice for a rough estimate of AF as a function oi q = |q| to take typical values for 
the other variables. We also neglect for now the imaginary part of AF . We assume < go < and, because the 
external electron lines usually represent members of a Cooper pair, we take |k|=fci? and fco = 0, measuring energies 
relative to Ep- Thus we have 



}_ [ [ dpxdpy 47re^ 
2^ J J (2^)2 (p-k)2+g^- 



AF/Fsardeen = 77— / / .^Ln/ 7Z , 2 ^^(P) P ' Qi^T^^L), (6) 



where k = fci?k and we have taken unit volume and set ?i = 1. I consists of six terms resulting from the po integrations 
which are displayed in the Appendix. For the p integrations we can set the Px axis along the direction of q and then 
the direction of k with respect to q is given by k^ = cos(9k)- For the present computation we have taken as typical 

values qn — O-lkp, kx = 0.8, uol — OAQEp, and uot — Q.MEp- 

For diagram Ic the complications due to the ionic corrections are absent and only two terms result from the po 
integrations. Then /S.Tcouiomb/^ Bardeen is given by Eq.(6) without the factor 1/20;^ and with I replaced by Icouiomb 
where 

l^Coulomb = —7 [ I ) 

£p+q - £p - qo - 10 Ep+g - Ep - qo + Id 

and f{£p) is the Fermi function at T = 0. The principal value integrations over px and Py must be done numerically. 

In Fig. 2 we show /S.T /T Bardeen and ^Tcouiombl^ Bardeen as fuuctions of q/kp for go = O.OSi^F (sohd curves) 
and go 0.02^^^? (dashed curve). There are two surprising results here. First, in 2D the vertex correction 
c ouiomb Bardeen IS already vcry large and, second, the effect of the ionic screening is drastic: One would naively 
expect the inclusion of the large ionic dielectric function to greatly reduce the correction. Except very near the sign 
change just below g — O.'ikp, this is not the case. A similar sign change was also found by Pietronero et al ||]. 

The structure in /ST couiomb/^ Bardeen deserves a brief comment. It arises from the complicated structure of the 
integrands and does not seem to have any obvious physical relevance. It is due to the first term in Eq. (7) which is 
singular at p = Ps(g, go) when go — Sp^+q ~ Sp^, i.e., when the external frequency can excite a particle-hole pair. The 
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integrand of the principal value integration in Eq. (6) then contains the singularity at p = Ps and a fixed (for k fixed) 
peak arising from the Coulomb factor. As q changes (for fixed k and qq), the singularity moves through the Coulomb 
peak causing alternately positive and negative contributions, producing the structure. This structure does not occur 
in AT/TBardeen becausc the integrand on the negative side of the singularity becomes positive a short distance from 
the singulatity and practically no negative contribution occurs. 



4) THE LARGE VALUE OF AT/TBardeen- 

The extremely large value of AT/TBardeen that we find here seems at first sight to be unphysical, since (to the best 
of our knowledge) it was not considered before. We suggest that it is physical. We are assuming here that our lowest 
order correction is not too strongly reduced by higher order vertex corrections, see discussion below. We believe that 
the reason why such an effect was not considered previously, is that in a homogeneous system, the dielectric constant 
e also enters into the expression for the bare vertex, which becomes: 

Ze^ I Q I /q^ £ , and not only into the expression for D,,^ ■ Our ansatz considers an extremely inhomogeneous system 

mm- 

While we consider effects due to phonons, we consider 3 entirely distinct phonon modes. Namely, the mode Q — qq 
which represents phonons associated with the momentum q of the bare vertex TbiocH, Fig. (la); and the modes cot 
and ul , which dominate the ionic dielectric constant. Since the modes are distinct, they involve motion of different 
atoms, in different regions of the unit cell. As mentioned previously, the mode ft could also be a spin fluctuation 
(or another type of Boson) . For now we assume it is a phonon involving motion of the planar oxygen; such motion 
could be longitudinal (along the Cu-0 bond) or transverse, in the a-b plane or in the c-direction. The frequency 
of the transverse motion is around 40 meV. This mode is seen as perhaps a McMillan-Rowell structure in tunneling 
experiments in Bi2Sr2CaCu2 08+5 |2^. The mode lot involves the motion of the Ba (or Sr) atoms in the c-direction. 
The mode ul involves motion of the apex oxygen, as well as motion of the planar oxygens. These modes are in the 
range 15-80 meV; lot is in the range 15-30 meV, and lol is in the range 40-80 meV. This distinction between the 
modes is important not only to avoid double counting, but also because it causes the shielding of the electron-electron 
interaction Dee by the ionic dielectric constant to be entirely different from the shielding of the electron-Jl phonon 
interaction (by the ionic dielectric constant, associated with the modes lot , i-^l)- We estimated this last shielding and 
found that it is small j^]. Therefore we do not introduce it into the present calculation. 

While the experimental determination of e{Lo) by the IR measurements |^,^ is definitive, it is instructive to consider 
possible microscopic causes for the anomalously large value of 

e . We believe that the large e may be related to a new degeneracy between the Zhang-Rice singlet and the anti 
Jahn- Teller triplet of the CuOs complex. Kamimura et al |25| calculated the splitting between the singlet and triplet 
as function of the occupation of the 3d shell of the chain copper, and found that for 0.55 holes there, these states are 
degenerate. This calculation is carried out using a quantum-chemistry algorithm. We can characterize this splitting 
by an effective Hubbard U, and write: Ue// = TSbare/s ■ TStare is the Hubbard U in a "normal" complex, i. e. several 
electron-volts. Thus the splitting of about one tenth of an eV calculated by Kamimura, indicates a value of e of about 
50. Also, this calculation shows clearly the local nature of this large dielectric constant. 

Anisimov et al |26j found that the singlet and triplet are close, by a rigorous LDA calculation. They also considered 
relaxation of the complex, i. e. different Cu-0 distances in the singlet and triplet states. They found a rather large 
relaxation, substantiating the ionic nature of the large dielectric constant. Stern et al found a near-degeneracy of 
the singlet and triplet states, near the Sr in LSCO by EXAFS measurements. Itai & Gatt carried out a quantum 
chemistry calculation, and showed the important role of the c-axis motion of the alkaline-earth atom. All these 
calculations show that the e's that describe the shielding of the electron-electron interaction Dee and the electron- 
planar oxygen interaction TbiocH are of an entirely different nature, and only the former one is greatly enhanced. 
Therefore the very large difference in their values is not " unphysical" . 

Since the diagram. Fig. le yields such a large result one can not rule out that further vertex corrections may also be 
important. In fact, the natural class of diagrams to include is the sum of "ladder" diagrams. Such diagrams represent 
nothing more than the spin susceptibility, or spin fluctuation propagator, and describe the (antiferromagnetic) spin 
correlations thought to be particularly important in the cuprate superconductors. Previous calculations p9| , ^ , which 
considered the effect of these contributions on the phonon frequencies and on the screening of the electron-phonon 
vertex in nearly antiferromagnetic systems, could be modified to include the ionic screening considered here. Thus 
the nearness of both a ferroelectric as well as an antiferromagnetic transition could perhaps be accounted for. 

Kim first considered effects of the proximity to a ferromagnetic or an antiferromagnetic transition, and showed 
that this can enhance the electron-phonon coupling constant. The enhancement arises both from softening of the 
phonon, and from an increase of the matrix element {k \ W \ k + q) . The first effect does not concern us here, (it 
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doesn't lead to an increase in in the strong-coupling case). The second effect bears some similarity to the effect 
that we consider here, except that we consider primarily the nearness to a ferroelectric transition. 

Kim considers the role of a sum-rule that restricts the magnitude of his effect. This sum-rule was also considered 
by McMillan ^ and originally proved by Heine, Nozieres, and Wilkins ||||. It states that 1/(0) = Z/N{Ep) . This 
rule applies to a homogeneous system. In an inhomogeneous system, it breaks down. Its origin is that when D^e is 
reduced from Dee by e, the bare vertex F is also reduced by e (Fig. lb) and these two effects cancel each other when 
q=0. (When g > 0, a large e actually diminishes V(q)). 

In an inhomogeneous system, the reduction of Dee and the reduction of F are given by two different dielectric 
constants. We showed pTIpP] that the e for F is close to one, while the e for Dee is large. This follows from the loeal 
nature of e 



We note in passing that the relevant values of q are on the order of 0.25A^ . pi| This is small enough that the use 
of the e derived for q=0 is justified. On the other hand, the screening length is about 4A which is small enough 
so the dielectric constant is local. 



5) BARDEEN'S VERTEX FOR AN lONICALLY SCREENED MEDIUM 

Bardeen calculated the sum over bubble diagrams (Fig. lb). The correction due to the "bubble" vertex is very 
large, and negative; i. e., for w = 0, the sum over the bubble diagrams reduces the vertex from the Bloch value 
'inZe^/q'^ by a factor of + DeeN{EF)] ~ + (qo/q)'^], as in Eq. (1). For q « kp, this factor is typically 
about 1/3, and for small q values, it is even smaller. This reduction accounts for the resistivity of monovalent metals 
being an order smaller than the value calculated assuming the Bloch value of the vertex. Now, to include the ionic 
screening, we should replace Dee m Fig. lb by Dee- Here we must set qo = in Eq. (2) to avoid double counting 
the electronic screening. Since Dee is so much smaller (at low frequencies) than Deej the vertex correction is much 
smaller, and when Dee can be neglected, the vertex becomes the Bloch vertex; i.e. it is considerably larger . Thus, a 
very large ionic dielectric constant restores the value of the vertex to the large original Bloch value at low frequencies 
(uj < LOT ~ 20meV). To illustrate this effect we consider T Bardeen replaced with T Bardeen where 



och J- Bloch 



Bardeen 



(8) 



In Fig. 3 we show T Bardeen/'^ Bloch and T Bardeen 1^ Bloch vs q/kp for e(0) = 40 and qo = /cf/2. It is clearly seen that, 
with increasing q, T Bardeen approaches T Bloch much faster then T Bardeen does. At the q values of interest here, the 
corrected vertex is quite close to the original Bloch vertex, e.g., for q = 0.2fci? we have TBardeen ~ ^-^^Bioch- 



6) CONCLUSIONS 

We have shown that the lowest order vertex corrections including ionic dielectric effects are large. In particular, we 
show that the electron- phonon vertex for small q and small lo is greatly enhanced over the "standard" Bardeen value. 
Indeed, surprisingly, the Coulomb vertex correction, l^Tcouiornb/^ Bardeem is already large. This may be related to 
the 2D nature of the system. We are, however, not aware of any explicit calculations of this quantity either in 2D or 
in 3D. It has apparently previously been considered phenomenologically to be included in the bare vertex. We have 
also shown that the inclusion of ionic effects also greatly affects the electronic screening. These corrections all tend 
to significantly increase the maximum superconducting Tc value. This may require us to reconsider the possibility of 
the high as being due, at least partially, to a phonon-mediated interaction. At the current level of approximation, 
however, one cannot make a definite statement in this respect. 

We point out that we have employed a " hybrid" formalism in which a 3D e is combined with 2d band structure and 
2D integrations. This is not done to make the calculations easier but to have a consistent picture in accordance with 
experiment. We do not believe that a pure 2D or a pure 3D senario is physically correct in the high-Tc superconductors. 

There are still a number of things to be done before a truly quantitative theory is obtained. We mentioned briefly 
in Section 4 that the large contribution of the lowest order vertex correction indicates that further vertex corrections 
may be important, in particular the particle-hole ladder diagrams that contribute to the spin susceptibility. At the 
least, one must consistently include such diagrams in the bubbles of the electronic screening diagrams as well as 
directly as vertex corrections. If the higher order vertex corrections can be summed to infinite order the result could 
turn out to be smaller than the large lowest order contribution. In any case, this difficult task must be carried out 
before quantitative comparison with experiment can be attempted. In addition, a better theory must go beyond the 
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present level of approximation in other respects, for example, the momentum dependence of TBardeen should also be 
included. More fundamentally, there is at present no microscopic, field theoretic derivation of the Lyddane Sachs 
Teller expression for e in Eq.(5), and in general it does not hold for a polyatomic unit cell. However, experimentally 
it applies extremely well for LSCO 0] and reasonably well for YBCO [||. 

We acknolowledge helpful discussions with J. Appel, A. Baratoff, J. Birman, O. Gunnarsson, V. Kresin, D. Manske, 
M. Peter, and K. Scharnberg. The work of M. W. is supported by the US-Israel BSF. 



APPENDIX: 

We display in detail the integrand of Eq. (6) . 

^ jcol - C.|)/(£p)/(£p+,) ^ K-^|,)[l-/(£p)][l-/(£p+,)] 

(wi - ep)(wL - Ep+q + qo) (wL + ep)(wL + Ep+q - go) 

^ (4-U;|)/(£p)[l-/(£p+,)] ^ (4-C.|)[l-/(£p)]/(£p+,) 

(£p - WL)(£p+g - Ep - go - iS) (ep + a;L)(ep+q - fp - go) 
_ [(gp+g - go)^ - - /(£p)]/(£p+g) _ [(£p+g - go)^ - ^\\!{ep)\\ - /(gp+q)] 

(gp+g - go - WL)(gp+q - gp - go) (gp+g - go + ^^L)(gp+g - gp - gO - iS) 

FIGURE CAPTIONS 



1. The various electron-Boson vertices discussed in the text: a) The bare (Bloch) vertex; b) The Bardeen vertex 
which includes electronic screening. Dee (q) = 47re^/g^; c) The lowest order Coulomb vertex correction; d) The 
electron- electron interaction including ionic dielectric effects (Eq.(2)). 

2. The vertex corrections ISX^ Couiomhl^ Bardeen and 1S£ jV Bardeen as fmictions of the external momentum q for 
external frequency go = 0.08i?F (solid curves) and go = 0.02i?F (dashed curves). 

3. The electronically screened vertex with Bardeen) and without (Psardeen) ionic screening. 
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